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Introduction 


Among  the  most  interesting  topics  in  chemistry  is  metastability  in  molecules.  Besides  telling 
much  about  the  chemical  bond,  it  offers  the  potential  for  identifying  molecules  that  can  have 
important  applications  from  high  energy  propellants  to  laser  .sources  [1].  Molecular  metastability 
may  be  achieved  in  two  predoininani  ways.  l‘irst,  we  can  speak  of  "geometric"  metastability. 
By  this  term  we  mean  some  orientation  of  the  atoms  in  a  molecule  that  is  generally  not  known 
experimentally,  but  that  has  stability  as  measured  by  it  being  a  local  minimum  on  the  potential 
energy  surface  (PES)  for  that  collection  of  atoms.  The  global  minimum  on  the  PES  will  be  lower, 
preferably  very  much  lower  if  we  hope  to  store  large  amounts  of  energy  into  the  metastable  form. 
Clearly  we  would  have  to  have  an  activation  barrier  sufficiently  high  that  separates  our  metastablc 
form  from  other  minima  and  particularly  the  global  minimum  so  that  we  can  trap  a  significant 
number  of  molecules  into  the  metastable  structure.  Such  a  generic  reaction  path  plot  is  shown 
in  Fig.  1, 


If  the  barrier  AEa  is  too  low  or  if  significant  tunneling  can  occur,  there  might  be  an  inadequate 
lifetime  for  the  metastable  state.  Alternatively,  if  there  are  lower  lying  states  that  offer 
decomposition  as  does  the  triplet  state  shown  in  Fig.  1,  then  spin-orbit  or  non-Born-Oppenheimer 
effects  can  cause  a  curve  crossing  again  jeopardizing  the  lifetime  of  the  metastable  state.  Thirdly, 
some  transients  can  rapidly  react  in  a  bimolecular  manner,  causing  a  short  lifetime  for  the  species. 
In  brief,  to  know  if  a  geometrically  metastable  molecule  might  be  synthesized,  we  need  much 
detailed  information  about  its  potential  energy  surfaces. 

'I’hc  .second  type  of  meiaslabiliiy  would  be  to  trap  a  molecule  in  an  electronic  excited  state. 
This  is  commonly  achieved  when  a  closed  shell,  singlet  molecule  is  excited  in  some  way  to 
a  triplet  excited  state,  e.g.  Such  a  state  is  forbidden  by  spin  .selection  rules  from  decaying  to 


2 


the  ground  state,  in  the  approximation  when  spin-orbit  and  other  spin  dependent  terms  are  not 
included  in  the  Hamiltonian  for  the  system.  This  leads  to  phosphorescence  and  in  the  absence 
of  collisions,  an  unusually  long  lifetime  for  the  excited  state  of  ~10'*  to  10s.  When  the  excited 
state  has  the  same  symmetry  as  the  ground  state  fluorescence  will  typically  occur  within  ~10'^s. 
Provided  we  can  form  a  geometrically  metastable  molecule,  we  can  always  attempt  to  further 
trap  the  molecule  into  an  excited  state  of  some  appreciable  lifetime,  making  the  two  ways  of 
achieving  metastability  additive. 

Since  metastable  molecules  cannot  be  easily  formed  and  studied  using  conventional  experi¬ 
mental  approaches,  predictive  quantum  chemistry  methods  have  an  essential  role  in  characterizing 
such  mctastable  species.  In  order  to  identify  high  energy  molecules  worthy  of  synthesis,  at  a  suf¬ 
ficient  level  of  approximation,  which  today  usually  requires  the  substantial  inclusion  of  electron 
correlation,  quantum  chemistry  can  provide  the  following: 


1.  Determination  of  critical  points  (minima  and  transition  states)  on  a  PES  using  analytical 
gradients  (first  derivatives). 

2.  Detailed  prediction  of  molecular  structures  (i.c.  bond  lengths  and  angles)  at  the  minima. 

3.  Evaluation  of  activation  barriers  at  the  transition  states  and  associated  lifetimes. 

4.  Unambiguous  characterization  of  the  critical  points  on  PES  as  minima,  maxima,  or  .saddle- 
points  (i.e.  transition  states)  by  evaluating  the  Hessian  (second  derivative)  matrix. 

5.  Prediction  of  heats  of  reaction  of  various  potential  metastable  forms. 

6.  Predictions  of  vibrational  spectra,  frequencies  and  intensities,  for  experimental  identification 
and  interpretation. 


7.  Determination  of  excited  state  energy  surfaces,  electronic  excitation  energies  for  UV/vis 
spectra  interpretation,  and  associated  lifetimes. 


8.  identification  of  other  electronic  states  that  can  affect  stability  of  mctastable  species. 


Combining  the  above  with  more  qualitative  theoretical  considerations,  such  as  what  decom¬ 
position  products  arc  allowed  by  symmetry;  and  what  kinds  of  bonding  patterns  can  be  expected 
to  result  in  viable  metastable  structures;  can  provide  a  microscope  that  allows  a  penetrating 
study  of  a  variety  of  possible  mctastable  molecular  forms  that  could  not  be  investigated  ex|>cr- 
imentally.  Furthermore,  we  can  investigate  situations  like  the  interactions  of  atoms  in  a  host 
lattice,  an  area  of  little  experimental  information.  All  our  work  is  made  possible  by  formulating 
and  implementing  state-of-the-art  ab  initio,  correlated  quantum  chemical  methods  based  upon 
coupled-cluster  (CC)  theory  [2-8]  and  its  finite-order  many-body  perturbation  theory  (MBPT) 
[3,4,9]  approximations.  Many  of  the  newest  methods  in  this  area  have  just  been  developed  in 
the  last  three  years  of  this  project. 

Largely  under  AFOSR  support,  we  have  been  developing  CCTMBPT  methods  since  the 
mid  seventies  [3,4].  These  methods  are  recognized  to  be  among  the  most  accurate  available 
in  quantum  chemistry  |31.  In  the  last  three  years  under  this  grant,  we  have  further  added 
to  the  theory  and  made  a  new  implementation  1 10]  into  a  highly  efficient  program  system  to 
provide  powerful  CC/MBPT  tools  for  antilytical  derivatives  relating  to  nuclear  displacement 
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(i.e.  gradient  and  hessian  matrices)  [11,12],  which  are  essential  to  locating  critical  points 
and  predicting  vibrational  spectra  including  that  for  open-shell  systems.  We  have  formulated 
and  succeeded  in  fully  incorporating  Abelian  symmetry  into  our  newest  CC/MBPT  methods 
[10].  All  of  these  recent  developments  have  been  implemented  into  a  general  purpose  computer 
program  (Advances  Concepts  in  Electronic  Structure  II,  i.e.  ACES  11)  to  evaluate  energies  and 
analytical  gradients  at  various  MBPT  [11]  and  CC  levels  [12,13].  The  timings  in  ACES  II  offer 
a  quantum  jump  in  efficiency  compared  to  other  widely  available  programs  like  GAUSSIAN  90 
that  provide  some  of  the  same  many-body  methods.  But  even  more  significantly,  many  of  the 
capabilities  of  ACES  II  are  unique.  Having  the  analytical  derivatives  themselves  is  only  pan 
of  the  problem,  however.  We  must  also  have  efficient  search  techniques  to  locate  minima  and 
especially  saddle  points  to  minimize  the  number  of  expensive  gradient  steps  required  for  systems 
with  several  degrees  of  freedom.  We  have  also  incorporated  several  such  approaches  into  ACES 
II,  in  hopes  of  making  searches,  particularly  transition  state  searches,  routine.  We  also  require 
second  derivatives  analytically,  or  as  numerical  differences  of  gradients,  to  characterize  critical 
points  as  minima  or  transition  states  [Il|. 

Furthermore,  we  include  in  this  report  .several  new  methodological  developments  which  will 
enable  us  to  better  describe  metastable  species. 


Metastable  Molecules 

To  date,  the  successes  of  our  effort  for  the  HEDM  program,  both  methodological  and  in  terms 
of  metastable  molecule  chemistry  have  been  many.  Considering  the  latter  first,  recognizing  the 
prospects  for  highly  energetic  metastable  species  based  upon  nitrogen,  we  initiated  our  study 
of  a  variety  of  such  potential  molecules.  These  included  N3H3,  which  is  isoelectronic  with  the 
well-known  species  ozone,  cyclopropane  and  propene;  tetrahedral  N4  which  is  isovalent  with  P,\ 
or  As4,  and  is  forbidden  by  symmetry  from  decaying  to  ground  state  N2  molecules;  octahedral 
Ng,  isoelectronic  with  cubane  and  similarly  forbidden  from  direct  decomposition  to  N2  ground 
state;  pcntazolc 


which  is  a  cyclic  aromatic  structure  i.sovaleiu  to  furan;  and  the  analogues  of  the  cyclopropenyl 
cation. 
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In  each  case  we  obtained  some  very  interesting  results. 

In  the  case  of  N3H3  [14]  we  obtained  the  structures  shown  in  Fig.  2,  that  consist  of  a 
triazene  form,  analogous  to  propene,  a  triimide  analogous  to  ozone,  and  a  completely  cyclic 
triaziridine  form  analogous  to  cyclopropane. 


Figure  2 


Each  of  these  was  found  to  contribute  a  metastable  minimum  on  the  global  N3H3  potential 
energy  surface  at  the  ab  initio  SCF  and  when  electron  correlation  (DZP-MBPT(2))  was  included. 
The  minima  were  fully  characterized  by  calculation  of  the  Hessian  (second  derivative  matrix) 
analyrically  at  the  MBPT(2)  level  [15],  which  also  provided  the  IR  spectra  of  these  unknown 
molecules.  At  this  level  of  approximation  the  vibnitional  frequencies  are  normally  within  5% 
of  experiment  and  the  intensities  typically  to  within  about  20%.  1  might  add,  the  character  of 
a  critical  point  on  an  energy  surface  can  change  depending  upon  whether  electron  correlation 
is  included.  E.g.  the  experimentally  observed  Cav  form  of  SiCa  [16]  is  not  a  minimum  at  the 
SCF  or  GVB  level,  but  is  at  the  MBPT(2)  level  [17],  so  this  is  a  critical  test  for  every  potential 
metastable  molecule. 

Although  sometimes  postulated  as  an  intermediate  in  inorganic  reactions,  there  is  no  experi¬ 
mental  evidence  for  the  existence  of  N3H3.  This  tends  to  be  a  recurring  theme  in  the  investigation 
of  metastable  molecules.  Yet,  the  accuracy  of  modem  ab  initio  correlated  predictions  is  suffi¬ 
ciently  reliable  that  the  molecule  should  be  capable  of  synthesis  at  least  in  matrix  isolation  and 
can  be  identified  from  our  predicted  frequencies.  Using  the  coupled  cluster  singles  and  doubles 
(CeSD)  plus  non-iterative  triples  T(CCSD)  model  [18]  gives  the  results  in  Tables  1  and  2.  Each 
N3H3  form  is  found  to  have  the  potential  for  storing  a  great  deal  of  energy.  Their  heats  of 
oxidation  correspond  to  specific  impulses  (Isp)  [19]  shown  in  Table  2  compared  to  456  for  liquid 
H2.  Furthermore,  the  formation  of  the  N3H3  molecule  is  pertinent  to  prospects  for  embedding 
N  atoms  into  an  H2  lattice. 


Table  1.  Relative  Energies  of  N3H3  Isomers  at  MBPT(2)  Optimum  Energies 

Relative  Energy  in  kcal 


Isomer 

MBPT(2) 

SDQ  - 
MBPT(4) 

SDTQ  - 
MBPT(4) 

CeSD 

CeSD  + 
T(CCSD) 

'Iriaziridine 

40.96 

41.21 

41.43 

40.84 

41.24* 

Triimide 

13.48 

2{).()1 

15.64 

20.58 

17.30 

Triazene 

0.00 

0.00 

0.00 

0.00 

0.00 

All  above  calculations  used  a  DZP  basis  set. 
.  *AH  =  87  kcal/mol  to  N2  +  H2. 

*AH  =  104  kcal/mol  to  N2  +  NH3. 


Table  2.  Specific  Impulse  of  Several  Molecular  Systems  (in  seconds) 


H2  +  O2 

456  (467)* 

N2H4  +  N2O4 

324 

N3H3  +  O2 

SCF 

CeSD  -h  T(CCSD) 

triazinc 

351 

347 

triimide 

382 

365 

6 


tiiaziridine 

399 

390 

N4(Td) 

520* 

489* 

N8(Oh) 

559* 

531  * 

♦The  formula  265\/A///Af  was  used. 


Our  next  effort  addressed  the  fascinating  question  of  whether  tetrahedral  N4  (tetraazatetra* 
hedrane)  constitutes  a  local  minima  on  the  N4  PES.  Symmetry  arguments  [20]  and  other  work 
[21,22]  suggested  as  much.  We  did,  indeed,  find  such  a  local  minimum  at  the  correlated  level, 
■characterized  it  via  correlated  second  derivatives,  and  predicted  its  vibrational  frequencies  [23]. 
Although  the  barrier  to  diss(x:iation  wa.s  dinicult  to  dctcnninc  definitively,  which  is  always  the 
case  for  a  Woodward-Hoffman  forbidden  process,  since  at  least  two  configurations  which  differ 
by  a  double  excitation  should  be  equally  important  in  the  transition  state  region.  Using  a  sym- 
metry  broken  multi-reference  CC  method,  we  obtained  a  good  estimate  that  the  barrier  was  more 
than  50  kcal/mol.  This  is  clearly  sufficient  to  make  the  molecule  experimentally  observable.  D. 
Yarkony  (private  communication)  also  tells  us  the  spin-forbidden  transition  to  a  triplet  surface 
has  a  barrier  of  at  least  30  kcal/mol  further  establishing  tetrahedral  N4  as  a  viable  metastable 
species.  In  establishing  its  existence,  however,  bimolecular  reactions  that  might  occur  cannot  be 
totally  disregarded.  We  also  located  four  other  minima  on  the  N4  surface,  including  a  C2v  and 
D2h  form  that  are  actually  even  more  energetic,  but  probably  would  not  have  as  high  a  barrier 
to  decomposition  (see  Fig.  3). 

Figure  3 

Relative  Energetics  of  Several  N4  Stationary  Points 
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Using  the  same  type  of  symmetry  argument,  the  Ng  (octaazacubane)  species  [22,24]  was 
similarly  identified  and  characterized  by  second  derivatives.  In  a  large  triple  zeta  double 
polarization  basis  (TZ2P)  the  energy  storage  at  the  CCSD  level  for  N4  and  Ng  is  respectively 
190.6  and  449.9  kcal/mol  resulting  in  Isp’s  of  489s  and  531s  (Table  2).  Both  exceed  the  H2  + 
O2  reference  value  of  467s.  Ng,  in  particular,  is  close  to  the  “revolutionary”  goal  of  6(X)s.  We 
also  looked  at  the  benzene  analog,  Ne,  in  this  paper  and  supported  it  not  being  a  local  minimum 
at  the  correlated  level. 


In  the  case  of  the  NBCH2  molecule,  six  local  minima  are  found  [25],  that  correspond  to 
the  variety  of  ways  the  atoms  can  be  attached  (see  Fig.  4).Although  not  as  energetic  as  N4  or 
Ng,  these  are  still  relatively  energetic  species.  Additional  pertinence  lies  in  their  relevance  to 
the  possible  local  structures  that  might  occur  when  embedding  N,  C  or  B  atoms  in  a  H2  lattice. 
Furthermore,  the  dimers  of  these  isomers  lead  to  natural  cyclic  structures  like 


Figure  5 


that  would  be  potentially  of  great  interest.  We  propose  to  study  such  dimer  structures  in  the 
next  year  of  this  research  program. 

The  above  nitrogen  ba.sed  species  are  energetic,  chemically  fa.scinating,  and  experimentally 
completely  unknown.  Consequently,  synthesis  in  suflicient  quantity  to  use  as  a  fuel  might  be 
difficult  to  achieve.  We  have  to  recognize  that  even  a  well-known  molecule  like  O3,  cannot  yet 
be  used  as  an  alternative  to  O2  in  an  actual  fuel  [26],  even  though  it  offers  a  greater  heat  of 
oxidation.  This  fact,  and  our  interest  in  providing  new,  but  realistic  high  energy  fuel  prospects, 
has  led  us  to  pentazole  based  molecules.  This  work,  which  is  unusually  promising,  has  just 
begun.  We  are  currently  making  a  theoretical  study  of  pentazole  derivatives  and  it  will  remain 
one  of  our  major  directions  for  the  next  year  of  this  program. 

Pentazole  (HN5)  has  the  fingerprints  of  a  metastable  molecule  since  even  though  it  offers 
a  perfectly  logical,  even  aromatic  electronic  structure,  it  is  unknown  experimentally.  However, 
unlike  N4  and  Ng,  derivatives  of  I IN5  like  phenyl  |)entazole  have  been  synthesized.  The  history  » 

of  this  molecule  is  intriguing.  Attempts  to  synthesize  it  go  back  to  1910  [27],  and  a  AgNs  is 
actually  reported  in  1915  by  Lisp.schiiz  128|.  Immediately  thereafter,  this  .synthesis  is  di.sputcd 
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in  the  strongest  possible  terms  by  Curtins,  et  al.  [29]  who  said,  “A  repetition  of  Lipschitz’s 
experiments  has  shown  that  all  his  observations  and  the  theoretical  conclusions  drawn  from  them 
are  erroneous  and  that  there  is  no  formation  of  pentazoles."  Several  further  attempted  syntheses 
were  not  successful  until  1957,  when  Huisgen  and  Uzi  wrote  a  paper  entitled,  “Pentazoles  I. 
The  Solution  of  the  Classical  Problem  of  Organic  Nitrogen  Chemistry”  [30].  This  paper  reports 
a  phenyl  pentazole  derivative.  Several  other  such  derivatives  have  been  made,  but  HN5  has 
never  been  synthesized. 

We  subjected  HN5  to  theoretical  analysis.  At  the  correlated  level  we  find  that  for  the  reaction, 
2HN5  — ♦  5N2  +  H2,  AE  is  -129.1  kcal/mol  resulting  in  an  Isp  of  346  as  a  monopropellant. 
Considering  that  derivatives  of  the  basic  structure  should  be  made  to  stabilize  the  pentazole  ring, 
we  can  consider  some  which  will  be  quite  energetic  (examples  include  nitropentazole. 


Figure  6 


and  a  bipentazole  analogous  to  biphenyl 


Figure  7 


and  possibly  a  dinitfamine  pentazole  (02N)N-N5,  and  numerous  metal,  MN5  structures,  where 
M  could  be  Li  e.g..  The  prospects  for  actually  making  and  using  pentazoles  is  promising. 
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Similarly  guided  by  realistic  synthetic  prospects,  we  have  also  considered  other  types  of 
high  energy  molecules  that  involve  Al,  B  and  Be.  The  primary  reason  for  investigating  sucli 
species  is  the  observation  that  the  heat  of  formation  of  molecules  containing  electron  deficient 
atoms  like  B  and  Be  is  exceptional.  In  fact,  Be(BH4)2,  shown  in  Fig.  8, 


Figure  8 


has  one  of  the  highest  heats  of  combustion  known,  -592  kcal/mol.  Although  known  for  50 
years,  this  molecule’s  structure  has  been  the  subject  of  controversy  throughout  that  period  [31 1. 
Some  of  the  structures  supported  by  different,  conflicting  experiments  are  shown  in  Fig.  9. 
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Figure  9 
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We  have  already  extensively  studied  Be(BH4)2  und  based  upon  making  accurate  predictions 
of  vibrational  spectra  with  correlated  methods  [32],  we  were  able  to  show  the  gas  phase  IR 
spectra  consisted  of  the  Dah  and  Dad  forms,  of  which  the  former  triply  bridged  form  solely 
remained  after  matrix  isolation.  There  is  still  a  discrepancy  with  the  apparent  observation  that 
the  molecule  has  a  dipole  moment,  since  neither  the  Dah  or  Dad  structure  can  have  one;  but  we 
attribute  this  observation  to  experimental  error.  Current  predictions  of  the  NMR  spectra  [33]  for 
the  Be(BH4)2  support  our  conclusions  of  the  structure. 

We  speculated  that  we  could  create  similar  molecules  with  exceptional  heats  of  combustion 
by  adding  a  third  BH4  group  to  the  basic  structure  of  beryllium  borohydride  with  a  Li  counter 
ion,  namely  LiBe(BH4)3.  This  molecule  is  unknown  experimentally.  In  preliminary  studies  we 
predicted  its  structure  und  find  that  the  lic(Ull4)7  anion  looks  like  a  propeller  (Fig.  10) 
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Figure  10 


and  finally  we  also  predict  its  vibrational  spectra  at  the  MBPT(2)  level.  Having  Be  in  a  propellant 
and  getting  complete  combustion  of  B  are  well  known  problems  to  the  explosives  community, 
but  the  exceptional  prospects  for  such  species  as  propellant  molecules,  perhaps  as  additives  in 
cryogenic  solids,  still  makes  their  investigation  particularly  pertinent  to  attempting  to  achieve 
revolutionary  Isp’s  in  realistically  synthesizable  structures. 

Another  example  of  a  readily  synthesizable  high  energy  material  is  Ai2B4H|8  which  has 
Dah  symmetry  (Fig.  11) 
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Besides  being  of  interest  because  of  its  high  energy,  a  correlated  description  of  this  molecule- 
using  a  DZP  basis  for  B  and  Al,  and  DZP  for  H,  requires  238  contracted  Gaussian  functions.  Not 
including  core  correlation,  there  are  224  active  orbitals.  We  determined  this  molecule’s  structure 
at  the  SCF  level,  which  is  as  shown  in  Fig.  1 1  above,  and  then  determined  its  energy  using  the 
CCSD  method  implemented  with  symmetry  in  ACES  II.  This  is  the  largest  CCSO  calculation  that 
has  been  reported.  The  calculation  requires  the  determination  of  two  million  cluster  amplitudes. 
Sustaining  an  execution  rate  of  290  megaflops  on  the  CRAY-YMP  the  CCSD  calculation  was 
accomplished  in  1303  seconds,  or  about  108s  per  iteration.  Including  all  steps  (integrals,  SCF, 
transformation,  sorting)  the  total  time  for  the  calculation  amounted  to  2973s  or  about  SO  minutes. 
The  AHf  was  determined  to  be  +65  kcal/mol.  From  the  CCSD  energies,  assuming  complete 
combustion,  the  AHrxn  -1600  kcal/mol  which  corresponds  to  an  Isp  =  528s. 

Methodological  Progress 

Several  recent  methodology  developments  have  been  accomplished  in  this  project.  These 
methods,  mostly  unique  to  our  effort,  make  it  possible  to  perform  high-level,  accurate  correlated 
calculations  on  much  larger  potential  metastable  species  than  was  previously  possible.  These 
include  the  following: 

1.  Direct  product  decomposition  approach  to  the  full  use  of  Abelian  symmetry  in  coUpled- 
cluster  and  MBFT  applications. 

2.  Restricted  open-shell  Hartree  Fock  (ROHF)  based  CC  and  MBPT  methods. 

3.  Analytical  first  derivatives  (i.e.  gradients)  for  open-shell  CC/MBPT  methods 
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In  the  study  of  metastable  species  high  symmetry  is  often  important,  since  for  such  examples 
symmetry  arguments  immediately  suggest  the  existence  of  barriers  to  decomposition.  Similarly, 
cluster  models  of  cryogenic  solids  will  frequently  exhibit  high  symmetry.  Most  ab  initio  quantum 
chemical  programs  that  can  be  applied  to  such  problems  do  not  take  advantage  of  point  group 
symmetry  once  electron  correlation  is  included.  We  developed  a  new  approach  to  readily 
incorporate  symmetry  for  sophisticated  methods,  termed  the  Direct  Product  Decomposition 
(DPD)  [10].  This  tool  underlies  all  of  our  new  program  system,  ACES  II.  To  cite  just  one 
example,  for  cubane  (CgHg)  a  6-3 IG**  basis,  “quadratic”  QCISD  calculations  (a  simplified 
CCSD)  which  is  the  highest  level  of  correlation  available  in  GAUSSIAN  90,  can  be  performed 
by  ACES  II  115  times  faster  than  by  GAUSSIAN  90,  and  the  whole  calculation,  SCF,  integrals, 
transformations,  etc.  40  times  faster  (see  Table  3).  The  time  required  in  our  calculation  is 
about  9  minutes  on  a  CRAY  YMP  compared  to  366  minutes  for  GAUSSIAN  90.  The  obvious 
impossibility  of  performing  such  calculations  for  large  molequles  and  basis  .sets  with  GAUSSIAN 
90  is  apparent. 


Table  3. 

Timings(s)  for  QCISD  calculation  for  cubanc,  6-31G'*‘  basis  set  (136 
basis  functions.  128  molecular  orbitals  included  in  QCISD) 


G90 

ACES  II 

RATIO 

Integrals  (VMOL“) 

799.8  5 

246.9  s 

3.2 

SCF 

1119.3 

26.0 

43.1 

lYansformation 

969.1 

45.6 

21.3 

Integral  Proc. 

16.0 

52.0 

0.3 

QCISD  Calculation 

19037.0 

168.2 

115.2 

QCISD  Iterations 

15 

11 

1.4 

QCISD  Cycle  Time 

1269.1 

15.3 

82.9 

TOTAL 

21941.2  s 

538.7  5 

40.7 

Almlbf  and  P.  Taylor 


Even  more  important  than  speed  arc  tlie  methods  available  to  treat  metastable  molecules. 
In  this  regard,  we  have  presented  the  first  generally  applicable  ROMF-CC  and  ROHF-MBPT 
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correlated  methods..  It  is  well  known  that  conventional  Hartree-Fock  (UHF)  reference  functions 
frequently  suffer  from  severe  spin  contamination.  The  alternative  of  using  a  restricted  open- 
shell  Hartree-Fock  (ROHF)  reference  largely  alleviates  this  problem,  but  the  more  complicated 
structures  of  the  ROHF  Fock  operators  makes  the  correlated  theory  difficult  to  formulate.  We 
found  a  simple  solution  by  formally  working  with  a  non-diagonal  spin-orbital  Fock  operator, 
where  all  new  terms  are  introduced  naturally  in  any  order  of  perturbation  theory  [34].  This 
differs  from  prior  work  of  Hubac  and  Carsky  [35].  in  that  we  use  the  Mpller-Plesset  unperturbed 
Hamiltonian  so  there  is  order-by-order  correspondence  with  the  usual  MBPT;  an  important 
element  when  open-shell  points  on  a  potential  energy  surface,  like  the  separated  atom  limit,  have 
to  be  compared  with  closed-shell  structures.  Our  approach  also  differs  from  that  of  Wolinski  and 
Pulay  [36],  in  that  ours  is  non-iterative.  This  makes  the  generalization  to  any  order  transparent. 

Finally,  our  ROHF-MBPT  approach  completely  eliminates  the  arbitrariness  inherent  in  spin- 
projected  MBPT(2)  results  [37,38].  The  latter  is  not  extensive,  as  required  for  a  many-body 
method.  Furthermore,  projected  teclmiques  do  not  eliminate  all  spin-contaminants,  as  in  our 
ROMF-MBPT  approach.  Lastly,  projected  MB1^'(2)  methods  do  not  have  natural  extensions  to 
higher  orders.  We,  however,  are  already  doing  ROIIF‘-MBFr(4)  and  analytical  gradients  through 
MBP1X3).  The  advantage  of  ROIIF-MBP'!'  over  other  techniques  are  readily  shown  in  Figure 
.  12,  where  the  electron  affinity  of  CN  is  plotted  versus  order  for  UHF-MBPT  and  ROHF-MBPT. 
The  latter  shows  a  very  smooth  convergence  to  the  infinite-order  CCSD  result. 


Figure  12 


Level  of  Theory 

Electron  alTuiily  of  CN  radical.  (O)  ROHF:  ( • )  UHF . 


Under  this  grant,  our  next  methodological  development  pertains  to  the  development  of 
analytical  gradients  for  all  levels  of  CC/MBPT  for  open-  and  closed-shell  systems.  We  presented 
the  theory  for  CC/MBPT  analytical  derivatives  some  years  ago  [39,40].  Now  we  have  augmented 
the  theory  and  implemented  it  into  ACIiS  II  to  apply  to  open-shell  systems  with  unrestricted 
Hartree-Fock  (UHF)  ]12],  re.stricted  open-shell  Hartree-Fock  (ROHF)  [41],  and  with  quasi- 
Hartrce-Fock  (QRHF)  references  ]13).  'I’hese  analytical  gradient  methods  apply  to  CCSD  [12], 
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CCSD(T)  [42],  QCISD  112J,  QCISD(T),  MBP1X2).  MBn’(3)  and  full  MBPT(4)  [43].  They  are 
the  only  open-shell  CC  analytical  gradients,  and  the  only  one  for  MBPT(4).  With  these  tools 
we  can  search  energy  surfaces  to  locate  the  extremum  points  and  transition  states  efficiently. 
Furthermore,  making  finite  displacements  of  analytically  computed  gradients  we  obtain  the  force 
constant  matrix  for  predicting  IR  spectra  and  for  documenting  the  existence  of  minima  and 
transition  states. 

Armed  with  these  developments,  we  are  in  a  unique  position  to  make  large  scale,  predictive 
ab  initio  studies  of  potential  metastable  molecules.  Our  state-of-the-art  methods  and  programs 
are  being  provided  to  other  AFOSR  investigators  for  their  use. 


i 
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